起伏地表弹性波传播的间断Galerkin有限元数值模拟方法   总被引：2，自引：1，他引：1       下载免费PDF全文

薛昭  董良国  李晓波  刘玉柱 《地球物理学报》2014,57(4):1209-1223

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.  相似文献   

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较     

董菲蕃  陈国兴  金丹丹 《地震学刊》2014,(2):154-160

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。  相似文献   

Application of perfectly matched layers in the transient analysis of dam–reservoir systems     

《Soil Dynamics and Earthquake Engineering》2014

The transient analysis of dam–reservoir systems by employing perfectly matched layers has been investigated. In previous studies, boundary conditions of the PML region in the reservoir have been neglected. In this paper, they are incorporated completely in the formulation. Moreover, a technique is introduced to involve the effect of incident waves caused by vertical ground motions at the reservoir bottom in the analysis. Performing several numerical experiments indicates that applying boundary conditions of the PML domain and utilizing the proposed method for vertical excitation cases reduce the computational cost significantly and make the PML method a very efficient approach for the transient analysis of dam–reservoir systems.  相似文献   

Detection thresholds in structural health monitoring     

《Soil Dynamics and Earthquake Engineering》2014

  相似文献   

Sampling efficiency in Monte Carlo based uncertainty propagation strategies: Application in seawater intrusion simulations     

《Advances in water resources》2014

The implementation of Monte Carlo simulations (MCSs) for the propagation of uncertainty in real-world seawater intrusion (SWI) numerical models often becomes computationally prohibitive due to the large number of deterministic solves needed to achieve an acceptable level of accuracy. Previous studies have mostly relied on parallelization and grid computing to decrease the computational time of MCSs. However, another approach which has received less attention in the literature is to decrease the number of deterministic simulations by using more efficient sampling strategies. Sampling efficiency is a measure of the optimality of a sampling strategy. A more efficient sampling strategy requires fewer simulations and less computational time to reach a certain level of accuracy. The efficiency of a sampling strategy is highly related to its space-filling characteristics.This paper illustrates that the use of optimized Latin hypercube sampling (OLHS) strategies instead of the widely employed simple random sampling (SRS) and Latin hypercube sampling (LHS) strategies, can significantly improve sampling efficiency and hence decrease the simulation time of MCSs. Nine OLHS strategies are evaluated including: improved Latin hypercube sampling (IHS); optimum Latin hypercube (OLH) sampling; genetic optimum Latin hypercube (GOLH) sampling; three sampling strategies based on the enhanced stochastic evolutionary (ESE) algorithm namely φ_p-ESE which employs the φ_p space-filling criterion, CLD-ESE which utilizes the centered L₂-discrepancy (CLD) space-filling criterion, and SLD-ESE which uses the star L₂-discrepancy (SLD) space-filling criterion; and three sampling strategies based on the simulated annealing (SA) algorithm namely φ_p-SA which employs the φ_p criterion, CLD-SA which uses the CLD criterion, and SLD-SA which utilizes the SLD criterion. The study applies SRS, LHS and the nine OLHS strategies to MCSs of two synthetic test cases of SWI. The two test cases are the Henry problem and a two-dimensional radial representation of SWI in a circular island. The comparison demonstrates that the CLD-ESE strategy is the most efficient among the evaluated strategies. This paper also demonstrates how the space-filling characteristics of different OLHS designs change with variations in the input arguments of their optimization algorithms.  相似文献   

How can the uncertainty in the natural inflow regime propagate into the assessment of water resource systems?     

《Advances in water resources》2014

The Canadian Rocky Mountain headwaters support the water resource systems of the Canadian Prairies. Significant variations in natural headwater contributions have been observed due to warming climate. Projecting future natural headwater flows under climate change effects, however, has large uncertainty. First, there are difficulties in climate modeling and downscaling in alpine regions. Second, streamflow modeling in mountainous areas is extremely challenging. There is therefore a need to understand the effects of uncertainty in the natural inflow regime, and in particular how this translates into uncertainty in representing the state and the outflow of water resource systems. Considering the Oldman River basin in Alberta, Canada, we synthesized different inflow regimes based on site/inter-site properties of the historical inflow regime. The water resources system was then conditioned on the synthesized inflow regimes to identify the mechanisms of error propagation from the headwater streamflows to the water allocations. The results show that the response of the water resource system to the uncertainty in the generated inflow regime depends on the system state, flow condition and the component of interest. Generally, the response of the reservoirs to the uncertainty in the estimated inflow regime is more significant in dry years, in particular during low flow conditions. The response at the system outlet is rather different, as the propagation of the headwater uncertainty is more significant during high flow conditions. Also, similar inflow estimates in terms of error and uncertainty may result in different error and uncertainty estimates in the simulated outflows; therefore, lower bias and uncertainty in estimating the regional inflow regime does not necessarily mean lower bias and uncertainty in simulating the streamflow at the outlet of the system. Our results provide improved understanding of uncertainty propagation through complex water resource systems, but also portray the need for better climate and hydrological modeling in the Rocky Mountains for improved water management in the Canadian Prairies, particularly in the face of uncertain climate futures. This will be crucial if the natural headwater inflows decline and/or the system faces drought conditions.  相似文献   

库木库里背斜活动性初步研究          下载免费PDF全文

孙浩越  何宏林  付碧宏  杨顺虎  白濱吉起  池田安隆  狩野谦一  越後智雄 《地震地质》2014,36(1):1-13

库木库里盆地位于青藏高原北缘,与柴达木盆地一山之隔,是二者的过渡地带,也是高原主体部分向NE扩展的前缘地区;现今构造表现为被3条大型活动构造带(走滑的阿尔金断裂带、东昆仑断裂带和逆冲的祁漫塔格褶皱逆冲系)所夹持。因此,该盆地对于研究青藏高原北缘的构造活动性、活动历史,探讨高原的扩展模式具有十分重要的意义。虽然库木库里盆地南、北两侧均发育活动性很强的大型走滑断裂,但是在盆地中央发育1条大型背斜,走向NWW-SEE,与祁漫塔格褶皱逆冲系和柴达木盆地内的褶皱构造走向一致,说明盆地目前遭受NNE向的挤压。通过对盆地地形横、纵剖面和阶地展布形态的分析,得出背斜有自西向东扩展变形的特征;野外调查和测年结果显示,背斜东段冰川融水形成了大型冰水扇,形成年龄为(87.09±2.31)~(102.4±3.7)ka,进而获得背斜东段自晚更新世以来平均隆升速率的最大值为(2.78±0.28)~(3.28±0.28)mm/a。库木库里盆地整体的活动性很强,在构造上与其北边的柴达木盆地类似,都受控于阿尔金断裂南侧的NNE向的区域挤压作用。  相似文献   

1D System identification of a 54‐story steel frame building by seismic interferometry          下载免费PDF全文

Mohammadtaghi Rahmani  Maria I. Todorovska 《地震工程与结构动力学》2014,43(4):627-640

A 54‐story steel, perimeter‐frame building in downtown Los Angeles, California, is identified by a wave method using records of the Northridge earthquake of 1994 (M_L = 6.4, R = 32 km). The building is represented as a layered shear beam and a torsional shaft, characterized by the corresponding velocities of vertically propagating waves through the structure. The previously introduced waveform inversion algorithm is applied, which fits in the least squares sense pulses in low‐pass filtered impulse response functions computed at different stories. This paper demonstrates that layered shear beam and torsional shaft models are valid for this building, within bands that include the first five modes of vibration for each of the North–South (NS), East–West (EW), and torsional responses (0–1.7 Hz for NS and EW, and 0–3.5 Hz for the torsional response). The observed pulse travel time from ground floor to penthouse level is τ ≈1.5 s for NS and EW and τ ≈ 0.9 s for the torsional responses. The identified equivalent uniform shear beam wave velocities are β_eq ≈ 140 m/s for NS and EW responses, and 260 m/s for torsion, and the apparent Q ≈ 25 for the NS and torsional, and ≈14 for the EW response. Across the layers, the wave velocity varied 90–170 m/s for the NS, 80–180 m/s for the EW, and 170–350 m/s for the torsional responses. The identification method is intended for use in structural health monitoring. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

Series expansion of complex ground amplifications with a sequence of simple transfer functions          下载免费PDF全文

Hiroyuki Goto 《地震工程与结构动力学》2014,43(13):1975-1997

Amplification of earthquake ground motions at actual deposit sites is an important factor to consider when assessing the risk of an earthquake disaster. In order to identify the amplification properties, several preprocessings such as the Fourier transform are required. I propose a series expansion of the amplification with simple ground transfer functions as a new preprocessing. I define a sequence of transfer functions based on a two‐layered structure excluding an internal damping and a function space spanned by the set of the functions. I mathematically prove that the function space is equal to L² space. This indicates that all the functions belonging to L² space, that is, an arbitrary ground amplification, have a unique series expansion. This expansion is applied to the physics‐based decomposition of the amplification. Some numerical examples indicate that the similarity between a target complex structure and a simple model is measured by the absolute value of each coefficient in the series expansion. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

地层错动引起的上覆饱和黏土层变形特性的离心试验研究     

骆冠勇  蔡奇鹏  吴宏伟 《岩土力学》2012,33(10):2985-2990

地震断层错动会引起上覆土层变形,从而造成断层附近的建筑结构、管线产生附加的变形和内力引起破坏。通过一个土工离心机试验分析上覆饱和黏土层在4步连续断层错动作用下的静力响应行为。着重分析断层错动引起的地层变形的范围、不均匀沉降区的分布特点、剪切裂缝在土层传播路径及地表开裂的位置等工程上重点关注的问题。得到以下几点认识：(1) 基岩断层错动引起的地层变形范围基本上不受基岩错动量大小的影响。(2) 断层错动引起地层的不均匀沉降区基本呈三角形分布,其地表宽度约为1倍左右的土层厚度。(3) 基岩错动引起的主剪切裂缝基本沿竖直方向向上传播,其传播距离取决于基岩错动量及土体的破坏应变。(4) 基岩断层错动在主剪区的下盘一侧边缘会产生张拉裂缝,且产生张拉裂缝所需基岩错动量远小于产生剪切裂缝所需的错动量。  相似文献